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By a canonical transformation of the three-band Hubbard model, we introduce an effective 
Hamiltonian for the propagation of two holes doped into the ground state of the Cu-0 plane. When 
the pair belongs to the B2 or 1 A2 Irreducible Representations of the d v Group, the bare holes do 
not interact by the on-site repulsion; the effective interaction between the dressed holes is obtained 
analytically in terms of renormalized matrix elements, and generalizes earlier findings from cluster 
calculations. The Fermi liquid is unstable and numerical estimates with reasonable parameters of 
the binding energy of the pair are in the range of tens of meV. Our scheme naturally lends itself 
to embody phonon-mediated and other interactions which cannot occur in the Hubbard model but 
may give important contributions. 



The discovery of high-Tc superconductivity §M has 
stimulated several serious and ingenious proposals for the 
pairing mechanism, ranging from enhanced BCS-type J2] 
and polaron ||] to a variety of purely electronic ones |J , 
just to cite a few of the less exotic proposals. These 
proposals have usually been regarded as alternative, but 
there is a clear need for a theoretical framework where 
they can coexist and can be compared on equal foot- 
ing, since all the serious mechanisms must convey part 
of the physics. Our first aim here is to propose such 
a scheme. Our approach started with cluster calcula- 
tions in the three-band Hubbard model; pioneers in this 
field explored the possibility that pairing could re- 
sult from repulsion; later it was pointed out JjJ that W=0 
pairs have intriguing properties. W=0 pairs are two-hole 
eigenstates of the kinetic energy Hq that are also eigen- 
states of the on-site repulsion term W with eigenvalue 0. 
The symmetry of the cluster is an essential ingredient, be- 
cause only fully symmetric clusters allow such solutions. 
The planar lattice structure is also essential, because no 
W=0 pairs occur in 3D or in a continuous model. The 
cluster calculations [p~0^p~2|| showed that W=0 pairs are 
the "bare" quasiparticles that, when "dressed" , become 
a bound state. That approach is inherently limited by 
the small size of solvable clusters, but allows a very ex- 
plicit display of paired hole properties, that even show su- 
perconducting flux-quantization Jl3| . The main indirect 
interaction through the background particles arises from 
the second-order contributions to the two-hole amplitude 
p3[ . Our second task here is to extend the theory to the 
full plane. The three-band Hubbard model Hamiltonian 
is H = Hq + W, and the independent hole hamiltonian 
reads, in the site representation 

Hp = y^£rfn.rf + y^gp^p + [c+c d + h.c] (1) 

Cu O n.n. 

where n.n. stands for nearest neighbors. The on-site 
repulsion Hamiltonian will be denoted by 
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W = yUin i+ ni-; 

i 



(2) 



where Ui = Ud for a Cu site, Ui = U p for an Oxygen. 
Omitting the band indices, we shall mean 

d[k] = \\k + , —k-\\ = cl + ct k _\vac > (3) 

to be a two-hole determinantal state derived from the k 
eigenfunctions. 

The point symmetry Group of the Cu-0 plane is C± v , 
and its character Table is shown in Table 1. We intro- 
duce the determinants Rd[k] = d[Rk],R e C± v , and the 
projected states 

* v [k] = ± £ x {v) (R)Rd[k} (4) 

where x^HR) i s the character of the operation R in 
the Irreducible Representation (Irrep) 77. In the non- 
degenerate Irreps, the operations that produce opposite 
Rk have the same character, and the corresponding pro- 
jections lead to singlets. Let Ri,i — 1, ..8 denote the op- 
erations of Ci v and fc, k! any two points in the Brillouin 
Zone (BZ). Consider any two-body operator O, which 
is symmetric (R^ORi — O), and the matrix with ele- 
ments Oi.j =< d[k] \R~l OR 3 \d[k'] >, where k and kl may 
be taken to be in the same or in different bands. This 
matrix is diagonal on the basis of symmetry projected 
states, with eigenvalues 

0{^k,k') = Y.^ ] ^ R ) R^ k ') ( 5 ) 

R 

where 

R (k,k') = (d[k]\6\Rd[k']y (6) 

Thus, omitting the k, k' arguments, we get in particular 

O (M 2 ) = O e + C2 + Ci + O c s 

-O ax - O ay - O a[ - O a , (7) 

O ( 1 S 2 ) = E + C2 - Ci - O c s 

-Oa, - O ay + O a[ + (8) 

If 6 is identified with W, since W E = W C2 = W ax = W ay 
and W C4 = W Ca3 = = W^, one finds W ( 1 A 2 ) = 
W ( 1 B-2) = . These are W=0 pairs, like those studied 
previously Jlq ]. One necessary condition for pairing in 
clusters is that the least bound holes form such a pair, 
and this dictates conditions on the occupation number. 
In the full plane, however, W=0 pairs exist at the Fermi 
level for any filling. 

Suppose the Cu-0 plane is in its ground state with 
Fermi energy Ep and a couple of extra holes are added. 
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In principle, by a canonical transformation one can ob- 
tain an effective Hamiltonian which describes the propa- 
gation of a pair of dressed holes, and includes all many- 
body effects. Let us see how this arises. 

The exact many-body ground state with two added 
holes may be expanded in terms of excitations over the 
vacuum (the non-interacting Fermi sphere) by a configu- 
ration interaction: 

|*o >= ^2a m \m > +Y^b a \a > + ^c^|/3 > +.... (9) 

m a 

here m runs over pair states, a over 4-body states (2 holes 
and 1 c-h pair), over 6-body ones (2 holes and 2 e-h 
pairs), and so on. To set up the Schrodinger equation, 
we consider the effects of the operators on the terms of 
|*o >• We write: 

H \m >= E m \m >,H \a >= E a \a >, ... (10) 
and since W can create or destroy up to 2 e-h pairs, 
W\m >= ^y m ', TO |m' > +J2\ a > W <x,™ 

m' a 

+ X>>W,j, ro . (11) 



V m i jTO vanishes for W=0 pairs in our model; however we 
keep it for generality, since it allows to introduce the ef- 
fect of phonons and any other indirect interaction that we 
are not considering. For clarity let us first write the equa- 
tions that include explicitly up to 6-body states; then we 
have 

W\a > = Y^\ m > + XV > Wa '> a 

m a' 

+Y)P > W , a (12) 


where scattering between 4-body states is allowed by 
the second term, and 

W\0 >= ]T K> W ™>,0 + l"> W ",f> 
m' a 

+J2\P'}W P ,, (13) 

0' 

In principle, the term can be eliminated by taking 

linear combinations of the complete set of states: when 
this is done, we get a self-energy correction to Ep and 
a renormalization of the vertices, without altering the 
structure of the equations. The Schrodinger equation 
yields equations for the coefficients a,b and c 

(E m — Eq) 

+^]am'K>,m' + ^baWm.a + ^Cf 3 W m .j3 = (14) 

m' a /3 
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(E a — E )b a 

+^a m 'W a . iTO ' + ^2b a ,W a . a > + ^cpW a .p = (15) 

m' a 1 (5 



(Ep - E ) cp + ^a m * A . rf + Yy^ W ^cc' = (16) 



where Eo is the ground state energy. Then, we exactly 
decouple the 6-body states by solving the equation for cp 
and substituting into the previous equations, getting: 



(E m — E ) a m + 

m' 



Vrr. 



E 



Eq — Ep 



W„ 



E 



Eq — Ep 



(17) 



(E a - E ) b a + J2a m ,W a 



+E & ° 



' E -Ep 



= 



(18) 



Thus we see that the role of 6-body states is just to renor- 
malize the interaction between 2-body and 4-body ones, 
and for the rest they may be forgotten about. If Eq is out- 
side the continuum of excitations, as we shall show below, 
the corrections are finite, and experience with clusters 
suggests that they are small. Had we included 8-body 
excitations, we could have eliminated them by solving 
the system for their coefficients and substituting, thus 
reducing to the above problem with further renormaliza- 
tions. In principle, the method applies to all the higher 
order interactions, and we can recast our problem as if 
only 2 and 4-body states existed. Again, the W a ^ a term 
can be eliminated by taking linear combinations of the 
a states: when this is done, we get a self-energy correc- 
tion to E a and a renormalization of the W m , a vertices. 
The equations become 



(E m — Eq) a m + Y^a m rV m . m i + ^b c 



(19) 



(E a - E ) b a + Y,a m ,W atTn , = (20) 

m' 

Solving for b a and substituting in the first equation we 
exactly decouple the 4-body states as well. The eigen- 
value problem is now 

(Eo — E m ) a m = 

5> m , {V m , m - + (m\S[Eo]\m')} , (21) 

m' 

where 



4 



/ lore 1 1 '\ ^ <m\W\a><a\W\m' : 
{m\S[E \ \m) = > . (22) 

a 

This is of the form of a Schrodinger equation with 
eigenvalue Eq for pairs with an effective interaction V+S. 
Then we interpret the wave function of the dressed 
pair, which is acted upon by an effective hamiltonian H. 
The change from the full many-body H to H is the canon- 
ical transformation we were looking for. However, the 
scattering operator S is of the form S = W e f f + F, where 
W e ff is the effective interaction between dressed holes, 
while F is a forward scattering operator, diagonal in the 
pair indices m ,m' which accounts for the self-energy 
corrections of the one-body propagators: it is evident 
from (21) that it just redefines the dispersion law E m , 
and, essentially, renormalizes the chemical potential, ft 
is important to realize that therefore F must be dropped. 
This happens also in the Cooper theory jlij , where an ef- 
fective interaction involving phonons is introduced via an 
(approximate) canonical transformation; the off-diagonal 
terms are kept, while the diagonal ones, representing elec- 
tron self-energy corrections, are dropped. Therefore the 
effective Schrodinger equation for the pair reads 

{H a + V + W eff )\a>=E \a> (23) 

and we are interested in the possibility that Eq = 
2Ep — A, with a positive binding energy A of the pair. 
The V interaction just adds to W e ff, and this feature 
allows to include in our model the effects of other pairing 
mechanisms, like off-site interactions, inter-planar cou- 
pling and phonons. 

We emphasized the fact that the canonical transforma- 
tion is exact because in this way our argument does not 
require U/t to be small. In practice, however, what we 
can do is to calculate W e f / neglecting 6-body and higher 
excitations, and keeping in mind that at least the struc- 
ture of the solution is exact when expressed in terms of 
renormalized matrix elements. So, in the following, we 
calculate the bare quantities. The a states are 3 hole-1 
electron determinants which carry no quasi-momcntum. 
We write 

\a >= | \\(k' + q + k2) + ,k 2 -,-q-,-k'_\\ > (24) 

where k 2 is the electron state and pedices refer to the 
spin direction; those with opposite spin indices contribute 
similarly and yield a factor of 2 at the end. From the 
interaction matrix element 

< \\(k' + q + k 2 ) + ,h-,-q-,-k'_\\ \W\d[s] >= 
5{q- s) U(q + k' + k 2 , -k', s, k 2 ) 
-6(k'-s)U(q + k' + k2,-q,s,k 2 ) (25) 

we find that the product in the numerator of (22) yields 
4 terms; two are proportional to S(p — s) and belong to 
F, while the cross terms yield identical contributions to 
W e ff. Using (5,6), we obtain the effective interaction 
between W=0 pairs: 
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< %[p]\W eff \%[s} >= 4 £ x M (R) E 



£/ (p, k 2 ,Rs + p + k 2 , —Rs) U (Rs + p + k 2 , — p, -Rs, fc 2 ) 
i (i?s + p + fc 2 ) + £ (s) + i (p) - e (fc 2 ) - E 

The sum is over occupied /c 2 with empty i?s + p + fc 2 . 
Note that W e jf does not depend on the sign of U. The 
diagonal elements recover the A expression derived from 
perturbation theory for clusters jl]| if Eq is replaced by 
the unperturbed eigenvalue 2e (p). 

Since actually V = in our model, we drop it in the 
following. The m and ml indices run over the projected 
eigenstates of the kinetic energy, and the k labels 

run over 1/8 of the BZ. We denote such a set of empty 
states e/8, and cast the result in the form of a (Cooper- 
like) Schrodinger equation 



2e (k) a (k) + E W ef f (k,k')a(k') = E Q a (k) (26) 



for a self-consistent calculation of Eo (since W e ff de- 
pends on the solution). Let Nc be the number of cells 
in the crystal. The U matrix elements scale as N^ 1 and 
therefore W e // scales in the same way. For an infinite 
system, Nc — * oo , this is a well defined, but quite in- 
tractable integral equation. The problem becomes dis- 
crete when working with a supercell of N$c x N$c = Nc 
cells, with periodic boundary conditions. However, A de- 
pends on U's and Nc in a complicated way, within the 
range of attainable supercell sizes. We define the Uniform 
Interaction Model (UIM) in which a constant effective 
interaction V prevails for k and k! in 1/8 of the empty 
part of the BZ. Having computed A for a given filling 
and Nsc according to (26) we can determine the value 
of V which gives the same A in the UIM. This will be 
V e f f , that is, the effective V of our theory, characterizing 
the strength of the attraction by a single quantity. Since 
V e ff is not strongly size dependent, its numerical conver- 
gence with increasing Nsc is observed. In addition, the 
UIM can be solved in the thermodynamic limit; writing 
E Q = 2E F - A we can estimate ^asympt for Nc — y oo . 
Indeed, the Cooper-like Schrodinger equation (26) with 
W eff = V > 0, leads to 



where p is the density of states, which is solved nu- 
merically. We use as input data the current estimates (in 
cV) t=1.3, e p =3.5, Ed—0, U p — 6s, Ud — 5.3s, where s is 
a scale factor induced by renormalization. 

At s=2.121, with E F =-1.3 eV, we get for 1 B 2 pairs 
the results shown in Table II. Here, n to t is the filling 
and V e ff is derived by comparison with UIM calcula- 
tions. We see that V e ff is fairly stable when increasing 



k 2 :e(Rs+p+k 2 )>E F 



e/8 




(27) 
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the supercell size and corresponds to A. asym pt values of 
about 20 meV. Convergence to the thermodynamic limit 
is achieved, and implies a Cooper-like instability of the 
Fermi liquid in this model of the Cu-0 plane. A asyrnp t 
values in the range of several tens to a few hundreds of 
meV are obtained with a scale s which is somewhat larger 
than 1; it could be that the current estimates of U's are 
a bit low; however, since the screening excitations arc 
explicitly accounted for in the Hamiltonian, it is reason- 
able that the input U's must be somewhat larger than the 
fully screened interaction. Moreover, contributions from 
phonons and other mechanisms can be included by a non 
zero V, and must be relevant for a comparison with ex- 
periment. Wc find that 1 A2 pairs are more tightly bound 
close to half filling, but x Bi pairs are favored when the 
filling increases. We get attraction and pairing at all fill- 
ings: since the present mechanism is driven by symmetry 
it works unless the system distorts. Further data will be 
presented elsewhere. 
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TABLE I. The Character Table of the C\ v Group 



N S c ntot A(meV) V eff (eV) A asympt (meV) 

18 1.13 121.9 7.8 41.6 

20 1.16 42.2 5. 9.0 

24 1.14 59.7 7. 28.9 

30 1.14 56. 5.7 13.2 

40 1.16 30.5 6.6 23.4 



TABLE II. Binding Energy of W=0 Pairs. 
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